Double-Exchange Model on Triangle Chain 
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We study ground state properties of the double-exchange model on triangle chain in the classical 
limit on t2g spins. The ground state is determined by a competition among the kinetic energy of 
the Cg electron, the antiferromagnetic exchange energy between the t2g spins, and frustration due 
to a geometric structure of the lattice. The phase diagrams are obtained numerically for two kinds 
of the models which differ only in the transfer integral being real or complex. The properties of the 
states are understood from the viewpoint of the spin-induced Peierls instability. The results suggest 
the existence of a chiral glass phase which is characterized by a local spin chirality and a continuous 
degeneracy. 

PACS numbers: 75.10.-b, 75.30.Et, 75.30.Mb, 11.15.-q 



I. INTRODUCTION 

The double-exchange mechanism H has received a spe- 
cial attention as one of the canonical mechanisms which 
explain the magnetism and transport properties of a class 
of transition metals. The Hamiltonian which describes it 
is known as the double-exchange model J, 3j . Assuming 
a two sub-lattice structure, de Gennes obtained the phase 
diagram as a function of temperature and hole concentra- 
tion, X Q. The ground state exhibits the antiferromag- 
netic (AF), ferromagnetic, and canted AF phases 
The disordered, two sub-lattice, and helical spin struc- 
tures and their degeneracies were discussed combining 
the two sub-lattice structure. These rich structures are 
induced from competition between the kinetic energy of 
the Cg electron and the direct AF coupling between t2g 
spins Him. 

On the other hand, studies without imposing the re- 
striction to the two sub-lattice structure have been re- 
ported. These results suggest the existence of the ground 
state without the translation invariance. An existence of 
the phase separation was reported mill mill. (For 
the model without the direct exchange between t2g spins 
see refs. m Q E E m ,18, 19, 20, 21].) The analy- 
sis based on the spin-induced Peierls instability ^3 
suggest the existence of a super-cell structure of the lo- 
calized spins. The spin state is distorted with a period 
which is commensurate with the Fermi momentum so as 
to open the Peierls gap, thus stabilizin g th e system. This 
is isomorphic to the Peierls instability |2J| . The only dif- 
ference is that the spin degree of freedom is distorted, 
instead of the lattice degree of freedom. 

In two or more dimensions, the degree of the phase fac- 



tor in the transfer integral plays an important role |25|, 
|2^ . For example, the generalized Peierls instability |23 
is expected to work and the ground state exhibits a quan- 
tum Hall effect 12^ if the state satisfies conditions psj . 
(See also ref |2^.) Actually the staggered 7r-flux state is 
stabilized at half-filling in the two dimensional square lat- 
tice 122, 113 . The anomalous Hall effect induced from the 
Berry phase was investigated The extensive 

investigations for the phase diagram were performed in 
two or three dimensions IH 0, |3l| and a flux state with 
super-cell structure was found away from half-filling. 

Recently the most intriguing is the study within the 
framework of Tsallis nonextensive statistics I^E IstI Is^ . 

In this paper, we study the double-exchange model on 
the triangle chain to investigate the effects of (i) the phase 
factor in the transfer integral and (ii) the frustration due 
to the geometric structure of the lattice. The triangle 
chain is constructed from the one-dimensional array of 
the triangles. (See Fig^) The triangle is the smallest 
unit for making a closed loop which is responsible for a 
non-vanishing flux up to the local gauge transformation. 
Then, the triangle chain has the simplest lattice struc- 
ture to study the effect of the phase factor. On the other 
hand, the geometric structure induces a frustration. The 
ground state is determined by the competition among the 
kinetic energy of the eg electron, the antiferromagnetic 
exchange energy of the t2g spins, and the frustration. We 
numerically obtain the phase diagram assuming the ab- 
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FIG. 1: Lattice structure of triangle chain. 
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sence of the phase separation. The ground state exhibits 
an insulator with a super-cell structure of the t2g spins 
for a wide range of the parameter space. The state has a 
continuous degeneracy which is a hallmark of the double- 
exchange model in one dimension. The spin structure is 
characterized by an incommensurate spin configuration 
[3^ and a local spin chirality. The spin chirality induces 
a non-vanishing flux to the Cg electrons. This is an in- 
teresting example of the finite spin chirality induced in 
electronic models without any assistance of other effects, 
such as the anisotropy in spin space, spin-orbit coup ling , 
lattice distorsion, and external magnetic field, etc.. [40j 



II. HAMILTONIAN 

The double-exchange model on the triangle chain is 
defined by 

L L/2 
i=l 1=1 

-\-J ^ ^ Si • Sj^ 



^■e- {{ig}} n ^ 4> and {{tf^jY n ^ 0, 

where A"^ is the complement of the set A. (See Fig|2Ia).) 
The latter, {{tfj}}" n {{ig}} ^ (f), is not obvious. (This 
holds even for the simple one dimensional chain with a 
periodic boundary condition in a finite system size. The 
effect can be attributed to the boundary condition by 
a local gauge transformation and becomes unphysical in 
the thermodynamic limit.) We show that there exist at 
least two examples in {{ig}}'^ n {{ifj }} for the model on 
a triangle lattice. This holds for any lattice except for 
tree. The first example is 



1.2 — — ''^1 — 



'-2,3 



3,1 



e 

= cos — . 
2 



(4) 



We set the relative angles, 8, small and the same for 
simplicity. (See Figl^b) for the corresponding t2g spin 
configuration.) We show that {{ig}} do not contain Q 
as a subset for any value of (0i, 0i, 92,4>2, d3,4'3)- Proof. 
We can set (^3,(^3) = (0,0) without loss of generality 



and get tW^ = 



cos(9/2) by setting 







(6,01, 0,4>2)- The rest parameter can be written by 



(1) ^1,2 



cos(e/2) cos(e/2) +e-^('f'^-'f'^') sin(e/2) sin(e/2), 



which never becomes cos(0/2) for ipi — (1)2 = 0. Because, 
in order to make 2 real, we have to set 4>i ~ (j>2 = nn 
where Ci^a is the annihilation operator of the Cg electron {ji =integer) which inevitably leads to if 2 = 1 or cos0. 



at the site i, and Si are localized {t2g) spins which are 
treated as classical vectors directed along {9i,(t)i) in the 
spherical coordinates. We denote the spin configuration 
by a set {9i, cpi} {i — 1,2, ■ ■ ■ , L). Moreover, J( > 0) is 
the direct exchange coupling strength between t2g spins. 
L (even) is the total number of the sites. The number of 
the triangle is L/2. In order to see the effect of the phase 
factor in the transfer integral, we study two kinds of the 
model. Their transfer integrals are given by 0, 0, 



This contradicts □ 

One might consider that this property is due to a finite 
Peierls phase along the path 1 — > 2 ^ 3 ^ 1, because 
the solid angle among Si, S2, and S'3 is quivalent to the 
magnitude of the flux penetrating the triangle. However, 
this naive observation fails. Even when the three spins 
are coplanar, there exist an element in {{^ij}}'^ n {{t^}}. 
The second example is 



,IJ.=R 



COS ^cos^+ e-*(^'-^.) sin ^ sin ^2) 



2 2 
6 



(3) 



for C- and R-models, respectively. Here Qij is the relative 
angle between the spins, Si and Sj. We impose the peri- 
odic boundary condition for the eg electron, c^+i = ci. 

For a given spin configuration, {6i, (jji} {i =1,2,- ■ • ,L), 
the transfer integrals, {^^l, are uniquely determined 
through or (Q. When we change the spin configu- 
ration, {9i,(j)i}, the range of {i^} forms a set {{i^}}- 
We study the difference of the ranges between {tfj} and 
{t^}, for the same domain, {{9i,(j)i}\ < 9i < 7r,0 < 
< 277, i = 1,2,- • ■ ,L}, or, in other words, the differ- 
ence between the sets {{i^^}} and {{^f^}}. For a system 
with two sites, the parameter space, 2 — zi,2, contains 
that of t{^2 = ki,2|: i-e- {^1^2} ^ {^1,2}- "^^^^ because 
*i;2 = ^1.2 = l-2i,2|e'*, $ = arg[tan(32:i,2/Kzi,2)], and 
we can set $ = without loss of generality. However, 
for any lattice ( except for tree ) the transfer integrals 0) 
and ^ have the mutually exclusive parameter ranges, 



ti2 — cos 



01, 



^2,3 = cos 



e 



2,3 



t^i — cos ■ 



e 



3,1 



, (5) 



-63,1 = 27r and if^ > 0- 



with the restrictions Oi,2 + 82,3 - 
(See Fig.|3c).) We show that {{i^j}} do not contain the 
set (O as a subset for any value of {9i, (pi, 02, 02, ^3, ^s)- 
Proof: We set (^3,03) = (0,0) without loss of generality 



{{tPi }}■ 



}} 





FIG. 2: (a) The relation between sets, {{tfj}} and {{t^^j}}. 
(b), (c) Examples of t2g spin configuration which realize the 
elements in {{tij}}'^ n {{^ij}}. (b) The eg electron gets com- 
plex transfer integrals by a motion 1-^2-^3^1 and (c) 
it gets real transfer integrals by the same motion. The three 
spins are coplanar in (c). 
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and get 

t^^ = 00362,3/2 = cos 6*2/2 > 0, 

t^^^ = 00363,1/2 = cos 0i/2 > 0, 

t'{2 = -cos 61,2/2 < 0. (6) 

This means that the one of the three transfer integrals 
in {tfj} should be negative. D The general proof is not 
easy. 



A. Several limits 

(1) For J = the spin configuration is expected to be 
ferromagnetic, because the double-exchange mechanism 
widens the energy band width of Cg electrons to gain the 
kinetic energy. The ground state is metallic. 

(2) For J = 00 or X = 0, 1, the spin configuration exhibits 
mutually 120-degree structure because of the frustration. 

(3) For X = 1/2, in this model the Wigner-Sites cell con- 
tains two sites and the state is a band insulator. 



III. METHOD 



B. Region JS^ 







To obtain the ground state we employ the method used 
in ref. [23|. When we assume an occurrence of the 2kF 
Peierls instability, the unit cell for the transfer integral 
contains 2q sites (or, in other words, q triangles) at the 
filling X = p/2q where p and q are mutually prime. The 
kinetic part in the Hamiltonian is reduced to 

L/2q 2q 
2q/2 

+ ^ '^2i-l,2i+lCN j+2i-l^N j+2i+l + 



We have two phases. 

(1) For X < 0.5, the state is the ferromagnetic metal 
(FM). The spin configuration is ferromagnetic and the 
amplitude of the transfer integrals becomes uniform. The 
unit cell of the transfer integral coincides with that of the 
Wigner-Sites cell. Then we have an empty upper band 
and a partially filled lower band. 

(2) For X > 0.5, the state exhibits an incommensurate |39| 
band insulator (BI-1). For x = p/2q, a ferromagnetic 2q- 
merization occurs. The spin configuration is completely 
ferromagnetic within the domain and incommensurate in 
long distance. All the energy bands are dispersive. The 
energy gap opens at the Fermi point. 



where 2q is the number of sites contained in the unit cell 
and L is assumed to be the multiple of 2q. To obtain the 
ground state, we numerically optimize the functional 



E{{9,,<j),})^Et{{9,,q^,}) 

2q 

+JS^ ^ [ cos 0i COS 6j + sin 0i sin 6j cos{(f>i — (pj 



where Et{{0i, 4>i}) is obtained by the numerical diagonal- 
ization of Q. 

In the calculation, the 2kp spin-induced Peierls insta- 
bility for the eg electron is assumed. For a sufficient large 
JS'^(^ 0.5) and x = 0.5,0.25, we confirmed that the as- 
sumption is valid in the following sense: We calculated 
the ground state energy for a more large super-cell struc- 
ture which is not characterized by 2kF modulation of the 
transfer integral and found that the energy is the same 
upto 10 digits as that of the 2kp, and the spacial period- 
icity is just the repetition of that of the 2kF. 



IV. PHASE DIAGRAM FOR C-MODEL 

We show the phase diagram for C-model in Fig|3fa). 
The phases are classified by the spin configuration, the 
modulation of the transfer integrals, and the position of 
the Fermi point. 



C. Region JS — * 00 

The spin configuration is nearly 120 degree structure 
and three spins in a triangle are coplanar. 

(1) For X < 0.5, the state exhibits the incommensurate 
band insulator (BI-2). All the sub-bands are dispersive. 
(This is in contrast to the simple one dimensional chain 
where all the sub-bands are dispersionless H^.) 

(2) For X > 0.5, the state is in the metallic phase (120M- 
1). The amplitude of the transfer integral is uniform as a 
unit of the Wigner-Sites cell and we have two dispersive 
energy bands. 



D. Region JS^ 



0.25 



The spin configuration has a rich structure because 
of the competition between the antiferromagnetic ex- 
change coupling and the kinetic energy of the eg elec- 
trons which favors the ferromagnetic spin configuration. 
The incommensurate insulator with dispersionless energy 
bands (BI-3) occupies a large region in the center of the 
diagram. Except for x = 1/2, the spin configuration and 
the transfer integral are incommensurate. (The property 
is the same that of the incommensurate gapful phase in 
the simple one dimensional chain j23j|.') In the masked 
region in the diagram, the three spins in each triangle 
are not coplanar. This phase is characterized by a finite 



4 



spin chirality and the non-vanishing flux when Cg elec- 
tron moves around one of the triangles. (The chirality 
is defined by ' ('^'2 x '5'3)}j, where j is the index 

of the triangle.) The state exhibits continuous infinite 
degeneracy. Therefore, the state is a chiral glass. Each 
degenerate ground state has a finite chirality. Because 
of the lattice structure, where the adjacent two triangles 
share only one site, the degenerate ground states with 
different magnitude of the chirality are continuously con- 
nected to each other. The continuous symmetry breaking 
is forbidden in one dimension. Therefore, in general, the 
chirality vanishes because the ground state is a superpo- 
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FIG. 3: Possible phase diagrams for C-model (a) and R-model 
(b) as a function of the electron concentration x and antifer- 
romagnetic direct exchange coupling JS^ between t2g spins. 
FM is the ferromagnetic metal phase. BFs are the incom- 
mensurate band insulators with incommensurate ferromag- 
netic polymerization (BI-1), with incommensurate spin (BI- 
2,4,5,6), and with dispersionless energy band (BI-3,7). 120M- 
1 and 2 are the metallic phases with nearly 120 degree struc- 
ture. AFM is the sub-lattice antiferromagnetic phase. The 
masked regions are the chiral glass phases. The calculation is 
performed at points where broken lines cross. 



sition of the states with different chiralities. 



V. PHASE DIAGRAM FOR R-MODEL 

The phase diagram for R-model is shown in Fig|3{b). 

A. Region JS^ ~ 

We have three phases. 

(1) ferromagnetic metal phase (FM) for x < 0.5, 

(2) incommensurate insulating phase (BI-4) for 0.5 < x 
< 0.85, and 

(3) metallic antiferromagnetic phase (AFM) for 0.85 < 

X. 

In BI-4 phase three spins in each triangle are coplanar 
within a numerical error. This property is different from 
that of BI-1 phase in C-model. The AFM phase is sub- 
lattice AF and the transfer integrals between the sub- 
lattices terminate. Then we have one cosine band and a 
decoupled dispersionless energy band. 

B. Region JS^ 00 

The spin configuration is nearly 120 degree structure 
and three spins in each triangle are coplanar. There are 
three phases: 

(1) metallic phase (120M-2) for 0.15 < a; < 0.35, 

(2) incommensurate insulator for < a; < 0.15 (BI-5) 
and 0.35 < a; < 0.85 (BI-6), and 

(3) incommensurate insulator with dispersionless energy 
band (BI-7) for 0.85 < x. 

The difference between 120M-1 and 120M-2 phases is 
as following: The transfer integrals are able to be nega- 
tive in 120M-1 whereas all the transfer integrals in 120M- 
2 are positive. The case holds between BI-5, 6, and BI-2 
phases. We could not find the difference between BI-5 
and 6 phases. In BI-7 the spin configuration and the 
transfer integral are incommensurate. In this phase all 
energy bands are dispersionless because some transfer in- 
tegrals vanish due to the spin configuration 

C. Region JS^ ~ 0.25 

The incommensurate insulator with dispersionless en- 
ergy band (BI-8) occupies large region in the center of 
the diagram except for x — 1/2. The transfer integral 
exhibits a modulation. The masked region is the chiral 
glass phase. (In the R-model, the transfer integral is re- 
placed by it's absolute value so the eg electron do not 
be affected by the phase factor.) In the vicinity of full 
filling, there are several phase boundaries of five phases, 
BI-4, 6, 7,8 and AFM. In low dope region, we have four 
phases, BI-5,6,8 and 120M-2. 



5 



-0.5 
-0.6 



\. y R-model 




\-. / y C-model 









0.2 0.4 0.6 0.8 1 

X 



FIG. 4: Estimates of the ground state energy for C-model 
(solid line) and R-model (broken line) as a function of electron 
concentration x at JS^ = 0.3. The error bars are smaller than 
the plotted points. 

VI. DISCUSSION 

We have obtained a possible phase diagram of the 
double-exchange model on the triangle chain at zero tem- 
perature. The ground state has a rich structure, which 
includes a band insulator with an incommensurate spin 
structure [s^ and a chiral glass phase. These states 
exhibit the continuous degeneracy. One of the origins 
of the degeneracy is due to the structure of the lattice 
where the adjacent triangles share only one site. The 
phase diagram is different from the simple one dimen- 
sional chain mainly in the following two points: (i) The 
structure of the phase diagrams of C- and R-models is 
different. (In the simple one dimensional chain, it is the 
same.) This is due to the enclosure relation between 
{{tg}} and {{tfj\}. The transfer integral of BI-2 be- 
longs to {{t%W n {{tf^Y, those of BI-7 and BI-6 (x > 
0.5) belong to {{tfj}Yn{{tfj}}, and those of FM and the 
states at x = 0.5 belong to {{tg}} D {{t^j}}. Reflecting 
this property, the different ground states are stabilized in 
R- and C-models in the same parameter in (x, J) space. 
For a sufficient large JS^, it is found that the ground 
state energy of C-model is lower than that of R-model 



for a; < 0.5 and vice versa for x > 0.5. The example 
is shown in Fig^ (Within our analysis, the energy de- 
pendence as a function of the electron concentration is 
concave near x ^ and 1. This suggests the absence of 
the phase separation. However, we need further investi- 
gation for the definite result.) Then we may be careful 
when we substitute or approximate ([^Jl by ijjji. (ii) For 
a small J, the global structure of the phase diagram is 
similar to that of the simple one dimensional chain. How- 
ever, when J is large enough, the 120-degree structure is 
dominant because of the frustration. Thus the AF phase 
which is appeared in the simple one dimensional chain 
vanishes. 

We use the numerical optimization method to obtain 
the ground state energy. In general, the method is known 
to be difficult to the function with many variables. We 
have two sites in the unit cell and two parameters, {9, </>}, 
for each site. Therefore, we have Aq independent parame- 
ters which should be optimized for an electron concentra- 
tion X = p/q. The inverse of the denominator is the reso- 
lution of the electron concentration in the phase diagram. 
Actually, for a; = 1/10 we have 40 parameters and the 
number of parameter is a realistic limit of our numerical 
calculation. We cannot examine the electron concentra- 
tion with more large denominator by the present method. 
The connectivity of the critical lines and the precise po- 
sition of the phase boundaries are not clear due to the 
above difficulty. The simple extension of the method for 
many parameters seems not to be realistic. 
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